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<N . 

We develop our novel model of cosmology based on the Bose-Einstein condensation. This 
model unifies the Dark Energy and the Dark Matter, and predicts multiple collapse of con- 
i-J^ ■ densation, followed by the final acceleration regime of cosmic expansion. We first explore the 

generality of this model, especially the constraints on the boson mass and condensation con- 
ditions. We further argue the robustness of this model over the wide range of parameters of 
^ , mass, self coupling constant and the condensation rate. Then the dynamics of BEC collapse 

and the preferred scale of the collapse are studied. Finally, we describe possible observa- 
^ ' tional tests of our model, especially, the periodicity of the collapses and the gravitational 

wave associated with them. 

§1. Introduction 

^ \ After WMAP, fundamental parameters in astrophysics, Hq, q, Qm, ^a, etc. 

^f) • have became more accurate and mutually consistent. Although the standard AC DM 

ly-^ , model can explain most of the observations, consistently with these parameters, it 

CD I still remains as a phenomenological model. In our previous paper, in order to im- 

' prove this model, we discussed a cosmological model in the framework of relativistic 

Bose Einstein Condensate (BEC) and gave a partial solution of "why now problem". 
There, the quantitative study has been limited within spatially uniform development 
^ ■ of condensates. In the present paper, we develop this cosmological model to include 

spatially inhomogeneous modes to describe the instability of BEC and to clarify the 
preferred scales of the collapsed objects. Furthermore, we establish the robustness 
of this model, especially showing a generality of the condensation strength. 

In general, BEC proceeds in bose gas of mass m and number density n, when the 
thermal de Broglie wave length A^s = \/27rh'^/ (mkT) exceeds the mean interparticle 
distance n^^^, and the wavepacket percolates in space, 

27r/i2n2/3 ^ ^ 

kT < . 1-1 

m 

On the other hand, cosmic evolution has the same temperature dependence since 
the matter dominant universe behaves, in an adiabatic process, as 



typeset using VTpT^.ds (Ver.0.9) 



2 



Fukuyama, Morikawa, Tatekawa 



pcxT^/^. (1-2) 

Hence if the boson temperature is equal to radiation temperature at z = 1000, 
for example, we have the critical temperature at present TcrUical = 0.0027K, since 
Tm oc and therefore, Ty/Tm oc a in an adiabatic evolution. Using the present 
energy density of the universe p = 9.44 x 10~^''g/cm'^, BEC takes place provided 
that the boson mass satisfies 

m < 1.87eV.*) (1-3) 

Conventional BEC is described in terms of the mean field which obeys the Gross- 
Pitaevskii (GP) equation 

ih^ = -^Ai; + V^ + g\^f^. (1-4) 

Here ijj {x, t) is the condensate mean field, and V (x) is the potential. The coupling 
strength g is related with the s-wave scattering length a as 

and therefore implies attractive interaction for a < 0. Originally, only positive 
value for a (and therefore g) has been considered for BEC, since negative value 
for a necessarily yields imaginary terms in the ground state energy and chemical 
potential. However, the appearance of the imaginary part in the energy simply 
implies that the ground state is unstable but BEC itself takes place, as a transient 
state, even in negative g.^^ This instability of BEC turns out to be crucial in the 
context of cosmology. 

The above GP equation is apparently non-relativistic. In the context of cosmol- 
ogy, we need a relativistic GP equation. The relativistic GP equation has a form of 
the Klein-Gordon equation with self-interaction and the Lagrangian is given by 

L = V^ (^g^''d^<P^d,(l> - m^4>^cp - ^(0^)2^ . (1-6) 

We discuss the metric tensor given by 

ds^ = (1 + 2^)dt^ -a^il- 2^)dx2 , (1-7) 

where = {t, x) represents the gravitational potential. The instability of this 
field has already been studied in^^ disregarding the cosmic expansion. We consider 
normal mass signature but the self-coupling is negative. The instability analysis is 
applicable for either signatures for and A. We will discuss about the difference 
of our formalism from the usual Higgs mechanism later. 

In section 2, the essence of the BEC cosmology is briefly summarized. Scalar dark 
matter plays an essential role in our scenario. We show in section 3 that BEC model 



*^ This constraint will be somewhat reduced later. See Eq.(32) 
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is quite robust, especially the final accelerated-expansion regime is always realized 
in wide range of parameters. In section 4, extending the phenomenological scenario 
reviewed in section 2, we argue the inhomogeneous modes and the instability of BEC 
based on the microscopic Lagrangian Eq.(l-6) with the metric Eq.(l-7). Section 5 is 
devoted to the miscellaneous observational problems to probe the remnants of this 
BEC cosmological model. 



§2. Basics of BEC cosmology 

We briefly describe the basic scenario of the BEC cosmology developed in.^)'^) 
The backbones of this model are 1. relativistic CP equation, 2. steady slow process 
of BEC, and 3. BEC instability which leads to the Dark Energy collapse. 

2.1. Relativistic Gross-Pitaevski equation 

From Eq.(l-6), the relativistic version of the GP equation*^ in the Minkowsky 
space, becomes 

A(P + m'^(f) + X{(p*(p)(f) = 0, (2-1) 



dt^ 

with the potential 



^ = mV</.+ ^(</.»2. (2-2) 



Substituting the decomposition of the classical mean field (j) = Ae^'^ and defining the 
momentum = —dp^S = (e, —p), where p = rwyv, 7 = (l — tP) , the relativistic 
GP equation reduces to the Euler equation for fluid: 

The energy-momentum tensor associated with Eq.(l-6) becomes 

T^. = = 2d^cl>*d,ct> - g^.{dcl)*d4> - m^ct>*ct> - ^(</.V)')- (2-4) 

For the isotropic relativistic fluid, it reduces to 

T^'' = dmg{p,p,p,p), (2-5) 
in the local rest frame. Here, the condensate part of p and p are given by 

p = = + n?(t>*(t> + ^{<t>*<t>f = 4>*^ + V, (2-6) 

and 

p = = = = - m^^v - -{(t>*(t>f = - v. (2-7) 



*^ We set c = 1, H = 1 hereafter. 
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2.2. Steady slow process of BEC 

We consider the cosmic evolution of various energy densities on average in our 
model, leaving aside the dynamics of inhomogeneous modes, derived from the micro- 
scopic Lagrangian (1-6) to section 4. Further, here we discuss the BEC cosmology 
phenomenologically, leaving aside the generality of the parameters to section 3. The 
whole evolution is given by the following set of equations.^) 

d \ ^ SvrG 

Pg = -3Hpg - Fpg, 

p^ = -6H (p0 -V) + Fpg - r'p^, 

PI = -3Hpi + r'p^. (2-8) 

Here F is the decay rate of the boson gas (i. e. uniform DM) into BEC, and i~" 
is the decay rate of BEC into collapsed BEC (i. e. localized DM)). The former F 
is a constant, but the latter F' appears only when the BEC satisfies the instability 
condition. These rates are transport coefficients which characterize the BEC phase 
transition although they should be, in principle, calculated from the Lagrangian 
Eq.(l-6) and the environmental informations. Here we fix these values phenomeno- 
logically in the present stage of our model*-'. As is easily seen from Eq.(2-8), DE(p^) 
and T>M.{pg + pi) arc intimately related with each other in our model; original uni- 
form DM{pg) condensates into DE(p^), and it collapses into localized DM(pi), which 
eventually becomes a dominant component in the total DM {pg + pi). 

There arc two relevant regimes of solutions to Eq.(2-8). One is (a) the over-hill 
regime, and the other is (b) the inflationary regime. The former appears when the 
condensation speed is high and the condensation overshoot the potential barrier, and 
the latter when it is low. A general evolution is a mixture of them; several regimes 
of (a) finally followed by the regime (b). Let us now examine these regimes. 

(a) Over-hill regime: This regime generally appears when the condensation 
strength is faster than the potential force and the condensation overshoot the po- 
tential barrier, especially in the earlier stage of the cosmic evolution. Actually, this 
regime is a fixed point of the first three equations in Eq.(2-8): 

(f> ^ CO, p^ ^ 0, Pg ^ 0, H ^ 0,a ^ a^, (2-9) 

The field goes over the hill of the potential, as in Fig. 1(a). 

The condensation speed F is fast at the first stage, and the bose gas density is 
simply reduced pg oc e"^* and the cosmic friction term becomes negligible. Then, 
Eq.(2-8c) yield 4> ~ V , and (p reaches singularity within a finite time. Since (p 
increases rapidly in the last stage of the fall, the BEC reduction rate —6H (p^ — V) oc 
Hcp^ dominates the BEC increase rate Fpg. Thus, we have eventually 0, and 

H ^ 0. However, actually, this virtual singularity is avoided by invoking the last 
equation in Eq.(2-8). 




*^ This is somewhat generahzed later in section 3.2. 
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Fig. 1. Schematic diagram of the evolution of BEC. (a) over- hill regime and (b) inflationary regime 
are depicted. The inflation appears as a result of the balance between the condensation and the 
potential force, V' = rpg/4>- 



(b) Inflationary regime: This regime appears when the condensation strength is 
weaker than the potential force, especially in the later stage of the cosmic evolution. 
This regime turns out to be a stable fixed point of Eq.(2-8): 

^ if*, 0,0-^0. (2-10) 

The BEC condensed field stops and stays at an intermediate position of the potential 
hill for ever, as in Fig. 1(b). This mechanism is a novel type of inflation, which is 
supported by the balance of the condensation speed (Fpg), and the potential force 

v = rpg. (2-11) 

Though both sides of Eq.(2Tl) exponentially decay to zero, the balance itself is 
automatically maintained*^ . 

In the actual universe, the above two kinds of regimes are realized successively. 
First, the over-hill regime repeats multiple times in general until the bose gas density 
is consumed and the condensation speed decreases. Eventually the condensation 
force balances with the potential force, and the final inflationary regime follows. 

In Fig. 2 (a), numerical results for the evolution of the cosmic energy densities are 
plotted. Black, dark gray, and light gray curves represent, respectively in this order, 
the cosmic energy densities of the bose gas (pg), BEC (p,^), and the localized objects 
(pi). In this example, four over-hill regimes are finally followed by the inflationary 
regime. 

Further, in Fig. 2(b), the evolution of the corresponding EOS parameter w = p/p 
is shown. It is clear that the BEC (p^) acquires the genuine DE property (i.e. 
w f» —1), only recently z < 3. For z > 3, p^p behaves as ideal gas (i.e. w 1). 
This is because the field 4> is in the stage of condensation and moving. Therefore 
it possesses kinetic energy and positive pressure. Thus, a genuine DE with w ~ — 1 
only appears for 2; > 3 despite we often call as DE and pg + pi as DM in this 

*^ This exponentially decreasing amplitude of the balance may lead to the instability of the 
inflationary regime and the autonomous termination of this regime, given some small external per- 
turbations. 
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Fig. 2. (a) Time evolution of various cosmic densities as a function of red shift z. This is a numerical 
solution of Eq.(2-8). The red line is Pg, the green line is p^, and the blue line is pi. Here we 
have set rh^ ~ 0.01, A = —0.1, F = 0.4. The variables with tilde are dimensionless ones defined 
in the text. Several BEC collapses take place, which finally followed by a phase with constant 
energy density (accelerating expansion), (b) Time evolutions of the w parameters. The green 
line represents = p^/ p^ as a. function of z. It is apparent that BEC behaves as ideal gas in 
the early stage and as cosmological constant in the later stage. Black solid line represents w of 
the whole system. 



paper. Thus, our model cannot be distinguished from the standard model with a 
cosmological constant, as far as we observe the cosmic evolution a{t). 



Now we explain the detail of our numerical calculations above and how to obtain 
physical scales from them. In Eq.(2-8) we simultaneously use parameters of quite dif- 
ferent orders of magnitudes, as a result of bridging microphysics to the macrophysics, 
such as the plank mass mpj, cosmic expansion rate H, and the boson mass m etc. 
First we normalize dimensional observables by two mass scales; m* for space-time, 
and mo for energy, whose orders are fixed later: 

T = m^t, H= — , p=— , m = — ,4)= — ,r = — , -T' = — , (2-12) 
m* rtiQ mo mo m=K m* 

where m^: and mo are related through the gravitational constant G or plank mass 
mpi : 

SttG 1 2 r ^ 

or mQ = mpim^, (2T3) 



"ic^ml mQ 



Then Eq.(2-8) is rewritten as 



= Pg + P4, + Pi 
Pg = -^Hpg - Fpg 



p/ = -3Hpi + r'p~^, (2-14) 
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where the prime means the derivative w.r.t. r. The ratio of the mass scales m^/rriQ 
appears only in the first term on the RHS of the third equation of Eq.(2-14), and this 
term docs not contribute since this factor is extremely small, m*/mo = mg/mp; <C 1. 

The scales rriQ and therefore m^, can be fixed, from our numerical calculation, as 
follows. The energy density of BEC is normalized as 

\ 2 X 

J + rh^^*^ + -{^*^f. (2-15) 

We have started our calculation, in the above example of Fig. (2), with ttt- = 0.1 
and obtained the numerical value for at present, p^o = 0.000015. We identify 
P<j>o = 0.73pcr0 with pcro = 9.44 x 10^^° g/cm^, from which it follows that 

mo 0.030eV, m = mmo = O.OOSOeV, and = 2.09 x 10~^^eV. (2-16) 

As we have mentioned after Eq.(2-8), DE{p^) and DM{pg + pi) are intimately 
related with each other in our model; DM(pg) condensates into DE(p0), and it col- 
lapses into localized DM{pi). Therefore, it is natural to expect that the amounts 
of DM and DE are almost the same. This is a kind of self-organized criticality^^ 
(SOC). Finally dominant DM{pi) is composed from the quantity which had been 
DE(p0) before. 

2.3. BEC collapse 

There is no degenerate pressure for bosons unlike fermions. Therefore, only the 
quantum pressure, expressed in the last term on the left-hand side of Eq.(2-3) , 
induced from the Heisenberg uncertainty principle can prevent the BEC to collapse. 
However, there is a maximum mass for this mechanism to work. The Compton 
wavelength of the object Xcompton = 27r/i/ (rnc) must be larger than the Schwarzschild 
black hole radius or its inner- most stable radius 3 (2GM/c^). This condition defines 
a characteristic mass scale 

^critical ~ fnli/m = Mkaup (2T7) 

only below which a stable configuration is possible for BEC. This structure is well 
known as a boson star.^) For example, the boson mass m = 10~^ eV gives McrUical = 
10^ ^Mq, which is almost the mass of Mercury. If the object mass exceeds this critical 
value, black holes are inevitably produced. There is no limit for the black hole mass. 
Thus, DE black holes of any size arc naturally produced in our model. However 
actually, shock waves would also be naturally produced in the collapse process of 
BEC. They convert the potential energy to heat and yields huge pressure, which 
may prevent the collapse to black holes. We will study the initial linear instability 
of BEC in the following sections. In any case, it will be true that many compact 
clumps (boson stars, black holes, hot clusters) are rapidly formed after the collapse 
of BEC. This fascinating scenario, early-formation of black holes and clumps, has 
been extensively discussed in^^ in other context. 
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§3. Robustness of the BEC Cosmology 

After a brief introduetion of BEC cosmology in the above, we would like to argue 
first the robustness and generality of the BEC model. In particular, we would like 
to clarify the condition under which the BEC phase is possible in the universe. This 
argument is deeply related with the basic parameters of the model: mass m of the 
boson and the condensation rate F. 

3.1. Mass constraints and the BEC condition 

Let us now consider how BEC is possible in the expanding Universe. In general, 
the charge density n of the boson gas is expressed as the sum of two contributions; 
the particle with positive signature and the anti-particle with negative signature: 

f '^'P^ r "'" ^1 

J (2^^e/3('^-/') - 1 ~ el^(^+^'') - r' ^^'^^ 

which is a function of temperature 1//3 and the chemical potential /Lt, and energy 
CO = y^p^+m? , in units of h = ks = c = 1. BEC takes place when /i = m and the 
critical temperature Tg or the critical density Uc is determined by setting so in the 
above. 

In the non-relativistic regime, i.e. m^, the above form yields, 

nc = C(3/2)(^y^' (3-2) 

or Tc = (27r/m) {n/( (3/2))^^^ . Below Tc or above Uc, the wave functions of individual 
particles begin to overlap with each other, i.e. the thermal de Broglie length exceeds 
the mean separation of particles, 

In this regime, the cosmic energy density of the non-relativistic matter has the same 
dependence on the temperature: 

/ 2^ \ 3/2 

n = no i^-J , (3-4) 

if we assume the entropy is conserved during the expansion and therefore p oc oc 
^13/2 'pijyg^ provided that the cosmic temperature had been once under the critical 
temperature at some moment in the non-relativistic regime, the universe would be 
always under the critical temperature and BEC can initiate all the time. 

The above fact sets the upper limit on the boson mass, for BEC to take place at 
the present universe. Let us suppose the boson gas was in thermal equilibrium with 
radiation in the far past. Suppose the boson had transformed from relativistic to non- 
relativistic at time with the scale factor a* = a (i*). The boson temperature Tb (t) 
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Fig. 3. Schematic diagram of the critical temperature Tcr and the cosmic temperature Tuniverse- 
In the non-relativistic regime (m < T), Tcr oc p^^^, Tuniv oc p^^^ . In the relativistic regime 
(m > T), Tcr oc p^^^ , Tuniv OC p^^*. The temperature evolution line for the universe is the same 
as the adiabatic expansion. 



at this time will be Tb (t*) = m = (t*). After that time, the boson temperature 
reduces as oc a~^, and the photon temperature T-y (t) as oc a~^. Therefore, for t > t*, 

Putting the present value of radiation temperature T^o = 2.73-fC into the above 
equation, we can estimate the present temperature of the boson gas Tb (to) = Tbo 
as 

Tbo = (^) r,o. (3-6) 

The present value of the critical temperature can be estimated from the present 
energy density po = 9.44 x 10~^^ g/cm^. The ratio of them is 

Tbo ^ Cm?"T^,om''\ ^3.^^ 
Tea 2ttpI'^ 

The requirement that this ratio is smaller than 1 yields the upper limit of the boson 
mass: 

m < ^'^f''';\ « 19eV. (3-8) 

It must be noted that this upper limit does not apply to the boson which has not 
been in thermal equilibrium with radiation in the past, or the boson composed from 
the pair of fermions. 

On the other hand in the ultra-relativistic regime, i.e. 3> m^, the critical 
density becomes 

nc = (3-9) 
and the cosmic energy density of the ultra-relativistic matter behaves as 

n = no(^) . (3-10) 
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Therefore, contrary to the nonrelativistic regime, even if the cosmic temperature had 
been once under the critical temperature at some moment in the ultra-relativistic 
regime, the boson temperature in the universe would eventually goes over the critical 
temperature and BEC would melt into thermal boson gas at that time. This trend 
is depicted in Fig. 3. 

3.2. Time dependent F and the robustness of the BEC model 

In the previous calculations,^-*'^) we assumed the condensation speed F is a 
constant parameter. However, this quantity F is not a simple term appearing in 
a Lagrangian, but a transport coefficient, which includes all the information of the 
many-body environment, e.g. temperature, density, fluctuations, etc. It should be 
calculated from the quantum field theory of finite temperature and density in the 
expanding universe, though such theory does not exist at present. Therefore we 
take the second best method, i.e. we try all possible time-dependent F. Although 
this does not specify F, we may establish some robustness of the BEC cosmological 
model. 

General transport coefficients would depend on the temperature and the density 
of the environment. Time dependence of such global parameters can be represented 
by a scale factor in the uniformly expanding universe. Therefore, possible time 
dependences of the parameter F will be exhausted by the inclusion of the scale 
factor, which normally depends in the form of power. Thus we assume 

F = ra{tf (3-11) 

where F and a are constants. Even in this case, the basic mechanism of the BEC 
cosmology does not change. Actually, the insertion of the expression Eq.(3-ll) in 
Eq.(2-8) is equivalent to recast the behavior of boson gas density, which is the source 
of condensation, as 

Pg<xa {t)~^ ^ pg(xa (3-12) 

while F F is still a constant. Intuition tells us that such change of the source gas 
density does not alter the scenario in essence. 

In order to check this intuition, we have performed several demonstrations in 
numerical methods. Results are in Fig. 4. 

In general, the reaction rate reduces when the temperature decreases. Therefore 
it will be natural to choose negative values for the parameter a. In all calculations 
with a = —1,-2,— 3, the qualitative behavior of the model, i.e. multiple BEC 
collapses followed by an inflation, does not change. Quantitative changes are the 
total number of BEC collapse and the identification of the present time. These 
results could be easily foreseen from the fact that the change Eq.(3-12) is equivalent 
to Eq.(3-ll). We have also performed positive values for a which may be less relevant 
in the actual universe. In this case, only a = 1 yields the qualitatively similar 
behavior, but the cases for a > 1 are not clear within our calculations probably 
due to numerical error. 
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(c) (d) 

Fig. 4. Time evolution of various cosmic densities as a function of red shift z. Same as Fig 2(a), 
but with time dependent condensation rate F — 0.1a(t)° where a = —1, —2, —3, and +1, 
respectively for (a), (b), (c), (d). Other parameters are set as ~ 0.01, A = —0.1. 

3.3. Numerical F and m -robustness of BEC 

Here we show the robustness of BEC on the relatively wide range of numerical 
values of condensation rate F and boson mass m, etc. 

The boson mass turns out to be 0.0024eV, 0.0030eV, 0.0041eV, 0.0033eV, 0.0019eV, 
0.0024eV, respectively for the parameters (a)-(f) in Fig. 5. The condensation rate F 
turns out to be 1.411 x W'^'^eY, 8.37 x W^^eY, 4.050 x lO'^'^eY, 1.0521 x lO-^igV, 
8.593 X lO^^^eV, 2.272 x lO^^^eV in the same order. The qualitative feature of our 
model did not change under the above variations of parameters. Further variation 
will be possible in principle, but our numerical code at present does not yield reliable 
results. This applies especially the late stage of the accelerating expansion, due to 
the exponentially reducing energy density. 

From the above results, we notice that the physical quantities m and F are not 
simply related with parameters m and F. Actually, the boson mass m only changes 
about factor 2 while the parameter mass rh changes about factor 16. On the other 
hand, the condensation rate F changes about factor 20 while the parameter F only 
changes about factor 4. 

Robustness of the boson mass value 0.003eV it O.OOleV may be somewhat inter- 
esting. However, we have not yet resolved the origin of this robustness; it may be due 
to the intrinsic nature of our model, or it may simply represent that our numerical 
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Fig. 5. The robustness of BEC on the wide range of numerical values of m, F, and A. (a) = 
0.01, f = 0.1, A = -0.1. (b) = 0.01, f = 0^4, A = -0.1. (c) = 0.01, f = 0.1, A = -1. (d) 
= 0.01, r = 0.4, A = -1. (e) = 0.04, T = 0.4, A = -0.1. (f) = 0.0025, T = 0.4, A = 
-0.1. 



range is too limited. Therefore, we will not specify the boson mass in the following 
arguments, but leave some possible rage of mass values open for further discussions. 

3.4. BEC instability and the collapse 

BEC collapse is a very complicated process. During the collapse of BEC, gravi- 
tational potential energy, GM'^/R, is released, where M and R are typical mass and 
scale of the collapsing region. If this collapse were free from shock waves, the col- 
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lapse were spherical symmetric and the collapsing BEC exceeded the critical mass, 
then most of the collapsing BEC would turn into black holes. However, such ideal 
conditions would never be fulfilled. 

If the BEC collapse takes place smoothly, and the process is adiabatic, then the 
temperature of BEC is given by 

T «i «i M^-TpT-\ (3-13) 

which is the same, if 7 = 5/3, as the behavior of the adiabatic universe in Fig. (3) , 
but the opposite time direction. The boson temperature is always below the critical 
temperature, T oc p"^^^, and BEC is maintained. 

If the BEC collapse takes place violently, and the gravitational energy always 
turns into uniform thermal energy, then the process is in virtual equilibrium, and 
the temperature of boson is given by 

T^^^^GMy^mp'/^ (3-14) 
R 

which b GMm/R — > m, in the limit that the system size approaches to 

the Schwarzschild radius R 2GM. This means that the boson becomes relativistic. 
Even in this case, the boson temperature is always below the critical temperature, 
and BEC is maintained. See Fig. 3. 

However in the real universe, shock waves are inevitably produced and fluctu- 
ations associated with the collapse would be enormous. As a result, some small 
portion of collapsed BEC terns into black holes and the rest portion would become 
normal boson gas. In any case, the universe becomes very clumpy at small scales. 

The collapsed BEC will gravitationally attract baryons to form a cluster, as in 
the standard CDM model. By contrast, at large scales, the potential is not affected 
and is the same as that in the standard ylCDM model. This point is further clarified 
below. 

We would like to restrict our study in the linear instability analysis in this paper. 
This is the subject in the next section. 

§4. Instability of BEC and large scale structure 

After examining robustness of the BEC cosmological model in the above, we will 
now argue how the instability of BEC, which is extending on the entire universe, can 
manifest to form localized structures. 

The BEC collapse inevitably takes place in the over-hill regime in our model. 
These collapsed components form localized objects and become direct seeds of the 
structures in the universe. Our special interest is the preferred scale of these struc- 
tures. Though the full dynamics of the BEC collapse would require involved nu- 
merical calculations, linear instability analysis is always tractable,, within semi- 
analytical calculations, to which we devote in this paper. 

The metric is chosen as 



ds^ = (1 + 2^)dt^ - a2(l - 2^)dx2 , 



(4-1) 
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where ^ = ^{t, x) represents the gravitational potential and a = a{t) the cosmic 
scale factor. The Lagrangian for the BEC condensate mean field cf) becomes, on this 
metric without a source term, 

L = a{tf{l - 2^) f (1 + 2^)-i|0|2 - a{t)-\l - 2^)-^(V0) • (V(/)t) 



~ a{tf{l - 4#)(^2 _ a(t)(V(^)^ - m'^a{tf{l - 2^)(^2 

-^aitfil - 2^),/.4 + (1 - 2<!>)aitfLg , (4-2) 

where the last line is the linearized form. The source term can be negligible in our 

analysis, since the BEC collapse takes place very rapidly compared to the slow steady 
condensation time scale. The equation of motion for the condensate becomes, up to 
the first order in 

+ 3^0 - ^V2<^ + m^(l + 20)(P + A(l + 2^)|(^|2<^ = . (4-3) 

The factor 1 + 2fp in the last term, which was not in,^) plays an essential role. The 
associated Poisson equation becomes 

= ATvGa^ + -(^ • -</.t + m^\4>f + ^|<^|^ - pol , (4-4) 
I a a 2 ) 

where po is the uniform background energy density. 

We now decompose the variables into the background component, with suffices 
0, and the linearly perturbed component, with suffices 1, as 

(f) = (f)o + (f)i , (4-5) 
# = + #1 . (4-6) 

The background is spatially uniform and only depends on time, while the perturba- 
tion is time and space-dependent. The background solution satisfies, from Eqs.(4-3), 
(4-4), 

^0 + 3^4 + (?ri^ + A|0oP)</>o = 0, (4-7) 
4>lh + (m^ + ^I'/'op) |</>oP - PO = . (4-8) 

In these Eqs.(4-7), (4-8), the variable 4>o can be assumed to be real without generality 
since all the coefficients are real. The equations of motion for the perturbations are 

= ^1 + 3-^1 - ^V^0i + m^{(f)i + 20i(j)o) 
a 

+ A (^<P14 + ^<Po<Pi + 2^1-^0 
V^^i = M^i + 4) + (m^ + A|0o|')'^o(<^i + 4)- (4-9) 
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The variable 0i is complex and <Pi is real. Therefore we replace them by three real 
functions: 

(t)i=x + iy, (4-10) 
^i = z, (4-11) 

where x = x{t,r), y = y{t,f), and z = z{t,f). These new functions x,y,z are 
decomposed into Fourier modes 

x = X exp (^nt + ik-fj , (4-12) 

y = yex.p (|j7i + ik ■ , (4-13) 

z = z exp (^f2t + ik ■ , (4-14) 

where x, y, z, f2, k are constants. Putting these into Eqs.(4-9), we have 

+ 3-n + ^ + m^ + 3X(pl J X + 2{n? + A^^)0o^ = , (4-15) 

l^f ^Z-n + %^m^ + M^y = Q, (4-16) 
\ a J 

(2(t)QQ + 2{rr? + \<t>l)(t>Qjx + ^;^z = Q. (4-17) 

Existence of a non-trivial solutions x,y,z requires that the above set of linear equa- 
tions are dependent with each other. Thus we have the condition 

n k2 \ k2 



a J AttG 



X r22^3-r2 + ^ + m2 + A(/>^ =0, (4-18) 
\ a / 

which determines the instability parameter as a function of the wave number k. 
If one of the solution f2 in this equation becomes positive for some k, such mode 
becomes unstable with the time scale f2^^*\ Since the shortest time scale for the 
structure of scale I = a/k to form is l/c from causality, such structure formation 
would be actually possible only if > I. Thus the structure of linear scale / is 
possible only if the condition k/a> f2 is satisfied. More precisely, the first structure 
formation takes place at the shortest time scale. This condition, by setting a as a 
positive constant smaller than unity, 

a- = f2 with < a < 1, (4-19) 

a 



*^ It should be remarked that even in the simple case a = ao = 1, (j)o =const, i?^ has positive 
real roots and is unstable for the case of < 0, A > unlike the statement in.''^ 
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uniquely determines the preferred linear scale o/fc* of the structure formed after the 
BEC collapse. 

The expression for ak/a = Q is solved for k, choosing the value associated with 
the most unstable mode among four solutions of f2 for Eq.(4-18), as 

k* I -^Iff + \/^tff + 647rG(l + a2) {-01^ - 2Km'^^^ + k2(/)6) 
~a ~ I 2(1 + a2) 

where adiabatic approximation H = 0, cpo = is utilized since the collapse time 
scale is much smaller than cosmic and condensation time scales. In Eq.(4-20), wc 
consider the regime so that wigyj = tu'^ — 3«;0q > 0, i.e. 0o smaller than the value at 
the inflection point of the potential V {(j)). This regime is first realized in the BEC 
condensation process. The above expression Eq.(4-20) can be further reduced to 

since in general, the present mass scales, i.e. of order eV, arc negligibly small in 
comparison with the Plank mass W'^^eV: nT-lfj ~ ~ fn^ ^ ~ G^^. 

A rough estimate of Eq.(4-21) and of the preferred scale is possible. Setting 
miff (^Q m^, Eq.(4-21) yields the linear scale Z* = a/fc* 

~ f 1 ~ io23cm « 30 f — V kpc. (4-22) 

\ m / m \m J 

The resultant mass scale associated with the BEC collapse, which took place at 
redshift z, would be 

M. = poz'^ll ^ 1.6 X 10^^ (^)"'mo, (4-23) 

in which strong mass dependence is apparent. Actually there are several different 
scenarios for the formation of localized structures, depending on the mass m of the 
boson. We suppose the first BEC collapses at about redshift z 20. 

(a) If the boson mass is about leV, then, the typical mass of the structure will 

be 

M* 1.6 X lO^^M© (4-24) 

which is of order of a galaxy. Some fraction of this mass turns into a black hole 
and the remaining boson becomes a hot gas surrounding the black hole. This is the 
typical structure expected from the BEC collapse in the present cosmological model 
for this boson mass. The detail of the mass of such black holes necessitates elaborate 
numerical calculations, on which we will report in a separate report in the future. 

(b) If the boson mass is about 10~^eV, then the typical scale well exceeds the 
size of the horizon, as easily observed from the power in Eq.(4-23). Thus, a structure 
is not formed in the early stage of the BEC condensation while the condition miff = 




k^ ^ Sm^^o 
a 



m, 



eff 



\ 
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m? — SncpQ > holds. In this case, BEC condensation further proceeds and crosses 
over the inflection point, beyond there ?7ie// becomes negative and we can no longer 
use the approximation ni^ff ~ ~ "^^^ Then, we have to go back to Eq.(4-20), 
which yields the solution 



Vl + q;2' 



(4-25) 



It means that the strong instability initiates immediately after the mean field crosses 
over the inflection point. Thus we define the time r as the elapsed time after crossing 
the inflection point. A structure of scale a/k is formed at around r = a/{ak). 
Since the BEC is very unstable and the time scale is short, we can expand ^ (r) = 
0inf + (^iniT + O (t^) , where (/)inf is the value of the condensation at the inflection 
point: i.e. m? = 3K0?jf. Utilizing the relations, (/c*/a)^ « \m? — ?iK(fP' \ /(I + o?) = 
6/t(/)infi;/'inf''"/(l + a^), we have = (6/t^inf0inf«^/(l + a^))^''^, and therefore k/a = 
(ar)"^ = (6K<^infi?!>inf/(Q!+a^))^''^- Putting approximate values = —6H {p^ — V)+ 
rpg Fpg, and pg ^ p,j,^ rri^cfy^, we have 

^ « (6«</>2r/(a + a3))V3 ^ (rn^rf^ « {{10-hV)ho-^\vf' « IQ-^^^V, 

" (4-26) 
which corresponds to the size U lO^km and M* Ri 10~^g. Therefore, BEC collapse 
cannot form macroscopic cosmological structures. These clumps would work as the 
ordinary DM and the scenario of large structure formation reduces to the standard 
model. 

Thus the above situations can be summarized as follows. We have two Jeans 
wave numbers: one is mP/nipi, which appears for the condensation 0o smaller than 
the inflection point, and (m^-T)^/^, which appears for <po larger than the inflection 
point. For the case (a) m leV, the former instability takes place and the latter 
has no chance to appears. For the case (b) m ^ 10~^eV, the former instability is 
not sufficient and the latter strong instability sets in to make the BEC collapse. 

(c) If the boson mass is far below 10^'^eV, then we can estimate the preferred 
mass scale utilizing the above scaling relation M oc m~^. For example, the boson of 
mass 10""^^eV would yield an object of the galaxy size, and lO^^^eV a cluster size. 
Because the mass is ultra-low, the boson would always be in the condensed phase. 
Therefore the BEC has a chance to form DM directly. This consideration naturally 
brings us to a popular idea that the DM around a galaxy or a cluster is formed 
from scalar field with ultra-low mass.^^~^^^ We will leave this possibility open in this 
paper, and proceed to the next subject; radiation of gravitational wave. 

§5. Observational remnants of the BEC cosmology 

We now turn our attention to possible observational remnants of the BEC cos- 
mological history, especially in the context of BEC collapses. Most prominent effect 
would be the emission of gravitational wave, which may be remaining as a fossil in 
our present universe. 
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5.1. Gravitational wave associated with the BEC decay 

In our model, universe repeats violent decay of BEC to localized objects in 
general. Associated with this process, gravitational wave is expected to be produced. 
The energy emission rate of the gravitational wave from the moving body can be 
calculated from the formula 



dE G fd^D, 



dt 45c5 V dt^ 

where 



(5-1) 



Dai3 = j P {^^aXfi - r'^Sais) dV (5-2) 



is the quadrupole of the whole mass distribution. Suppose the object of linear size 
R collapses with the typical speed v. Then the total energy emitted would be the 
integration of the above formula during the collapsing time scale. 

This is roughly the gravitational potential energy of the extended object multiplied 
by the 'efficiency' {v/cf' . We now estimate possible remnant gravitational wave in 
the present background sky, for several cases bellow. 

(a) If we adopt the mass of the boson is about leV, then from Eq.(4-24), the 
preferred scale is i? = 10^^ cm, M = 1.6x lO^^M©. If we tentatively assume v = c/10, 
then 

E « 10^^ erg (5-4) 

which should be compared with the total rest energy of a star: Mq = lO^^'erg. If we 
further assume that the first BEC collapse took place at redshift z 20, then the 
present energy density of the gravitational wave becomes 

P9r,z=0 = Pgr,z=20 ( 20)-' = ^ (20)-' = lO^^l-^ (5-5) 

R-^ cm-^ 

which should be compared with the total energy density at present, per = 10~^^g/cm^ : 
10"*^ erg/cm^, thus 

^9r,z=0 = 10-'^ (5-6) 

Since the strain h associated with the gravitational wave is related with 

where lo is the frequency of the wave, we have 

h ^ 10~^^ for u « (30kpc)~^ = lO^^^Hz, (5-8) 
h 10-^^ for u lO^Hz, (5-9) 
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where the frequency is estimated by naive extrapolation. This should be compared 
with the present limit of the gravitational wave h 10"^^ for uj k. lO^Hz. 

The gravitational background formed during the inflation is estimated^^^ as 

Qg^ = {H/mpif Qr ~ (10^^)^ 10^^ = 10"^^, (5-10) 
and the associated strain is 

h ^ 10"^^ for Lo ^ lO^Hz. (5-11) 

The gravitational background formed during the oscillation of the cosmic string^^^ 
is 

f2g^ = 100 (G/u/c^) ~ 100 • 10"^ • 10""^ = 10~^. (5-12) 

(b) If we adopt the mass of the boson is about lO^'^eV, then from Eq.(4-26), the 
preferred scale is i? = lO^cm, M = 10~^g. Then the present energy density of the 
gravitational wave becomes 

Pgr,z=0 = 10-51 ^ , (2gr,,=o = lO'^^ ^ (5. ^3) 

cm'^ 

mainly at frequency oj ~ (lO^cm) ^ ~ lO^Hz, which is totally small and would be 
never be detected. 

The overall amount of energy density, as estimated in Eq.(5-6), will not affect 
the present standard cosmology. However, the strain, as estimated in Eq.(5-8), may 
have a chance to be detected as well as the case of gravitational wave formed in the 

inflationary stage. 

5.2. log-z periodicity 

The BEC collapses do not take place randomly but they are periodical events in 
the logarithm of cosmic time. As argued in section 2, the bose gas density is simply 
reduced as pg oc e~^* in the over-hill regime since the condensation speed F is faster 
than the cosmic dilution time scale. Therefore pg is simply transformed into p^. Just 
after each collapse, new BE-condensation always begins from = until it reaches 
some critical value p'^ = 0(1) Vmax ~ m'^/ (— A). Therefore the condensation energy 
density behaves [pg (to) — Pg (t)] modulo p"^ and 

P<t> (t) ~ [pg (to) - Pg (Olmod ~ [p9 (^o) (l - e"-^*)] ^^^^ , (5-14) 

where to is the time when the first condensation begins. 

Accordingly, we expect that each BEC collapse takes place after the time interval 
At from the preceding collapse at time t, where At is determined by the condition 

P^ = Pg (to) (e-^* - e-^(*+^*)) (5-15) 

This implies that the occurrence of the BEC collapse is periodic in the logarithm of 
time, log(t). If the cosmic expansion is power law in time, i.e. a{t) oc t^onst^ then 
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log(i)- periodicity directly implies log(a) and log{z) -periodicities. For example, 
in the typical numerical calculation in Fig. (2), BEC collapse takes place at 2; = 
{33.2, 20.2, 13.3, 8.89, 5.50}, which is almost log-periodic. 

This log(z)-periodicity is a general consequence of our model, provided that the 
over-hill regime and BEC repeats several times. Furthermore, this periodic BEC 
collapse may leave its trace in the non-linear regime in forms such as the discrete 
scale invariance or the hierarchical structure in the universe, provided the scale is 
appropriately chosen. 

As was discussed in the previous section, the preferred scale of BEC collapse 
only depends on the basic parameters of the model; /* ~ mpi/m? for m > leV and 
h ~ (m2r)V3 for m < lO'^eV. Then the sequence of the cluster produced through 
BEC collapses has the series of mass proportional to z"^, where z is the redshift of the 
collapse. Thus, the cluster mass also has the log-periodicity, and the largest cluster 
is formed at the first BEC collapse. The most interesting case would be when the 
boson mass is about leV. Then the hierarchy of galaxies are formed, for example 
in Fig.(2), M/Mq = {1.6 x 10^\ 3.6 x 10^°, 1.0 x 10^°, 3.1 x 10*^,7.2 x 10^}. Detail 
of the arguments on the observational size is a future problem. If the boson mass is 
about 10~^eV, then the hierarchy would yield no interesting scales with respect to 
the large scale structure formation. 

§6. Summary 

We have developed the cosmological model based on the Bosc-Einstein Conden- 
sation (BEC) from various points of view. This BEC cosmology is characterized by 
(1) the unification of DE and DM, (2) their mutual conversion, (3) quantum me- 
chanical condensation as a novel phase of DE, (4) violent collapse of DE, (5) log-z 
periodicity of the DE collapses, (6) black hole formation from DE, (7) formation 
of localized objects in high redshift regime, (8) inevitable final phase of accelerated 
expansion, etc. These have been briefly explained in section 2. We have examined 
this model in detail especially with respect to its robustness and instability in this 
paper. 

We have first examined the robustness of the model in section 3. By assuming 
thermal equilibrium of the boson field and the radiation in the early universe, we set 
the upper limit of the mass for the boson, which turned out to be about eV. Thus we 
have showed that BEC takes place very naturally in the universe. Moreover, if the 
boson has not been in equilibrium with radiation in the past, then even any value of 
mass is allowed. Next we demonstrated that the BEC takes place in the wide range 
of microscopic parameters of boson mass m, self coupling A, and the condensation 
rate F. Furthermore, we have revealed that the time dependence of F does not 
qualitatively affect the BEC model. Thus we have been able to show the robustness 
and the naturalness of the cosmological BEC model. 

Then we have examined the instability of BEC in section 4. We have calculated 
a preferred scale of the structure formed after the BEC collapse associated with this 
instability. This scale turns out to be quite sensitive to the mass of the boson. If 
the boson mass is about leV, the preferred scale is Z* ^ mpi/m^ and it is about a 
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galaxy size. If the boson mass is about 10~^eV, the scale is U k, {w?rY^^ and it 
is about a gram. If the boson mass is far much smaller, there is a possibility that 
DM is formed as BEC, and the preferred scale can be of galaxy or cluster size. We 
have also estimated the amount of the remnant gravitational wave associated with 
the BEC collapse. It turns out to be marginally observable in the near future if the 
parameters of our model is most optimized, otherwise it is simply too small. 

The present results as a whole, suggest that the boson mass is probably of order 
to leV. These mass scales are so tiny that no empirical evidence for such 
boson has been found yet. However, it is potentially interesting that these mass scales 
are the same order as the neutrino masses. Therefore it would be natural to consider 
that the boson particle is a composite of neutrino-neutrino(or neutrino antineutrino) 
pair^^)'^^) though this requires further investigations including the problem of how 
Fermi surface can be stable in such tiny mass. We would like to further extend the 
cosmological BEC model in our next paper. 
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